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Introduction and Objective 

Traj Features: 

• Program used to simulate atmospheric flight trajectories of entry capsules [1] 

• Includes models of atmospheres of different planetary destinations - Earth, Mars, Venus, Jupiter, Saturn, Uranus, Titan, ... 

• Solves 3-degrees of freedom (3DoF) equations for a single body treated as a point mass 

• Also supports 6-DoF trajectory simulation and Monte Carlo analyses 

• Uses Fehlberg-Runge-Kutta (4 th -5 th order) time integration with automatic step size control 

• Includes rotating spheroidal planet with gravitational field having a J 2 harmonic 

• Includes a variety of engineering aerodynamic and heat flux models 

• Capable of specifying events - heatshield jettison, parachute deployment, etc. - at predefined altitudes or Mach number 
Has material thermal response models of typical aerospace materials integrated 
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Modify trajectory simulation tool, Traj, to make it suitable for meteor entries including mass loss & fragmentation 
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Modifications Made to Traj for Meteor Simulation 


NASA's Galileo probe to Jupiter only one that experienced significant mass loss 

Entry capsule was a 45° sphere-cone with fully-dense carbon phenolic as heatshield material 

M. Tauber et ol. [2] developed JAE code for simulation of Galileo probe (Jupiter entry) 

JAE logic incorporated into Traj 

— Sphere-cone shape replaced by sphere 

— Mass loss equation of meteor physics used 

— Allow input specification of heat of ablation, Q 

— Allow heat transfer coefficient to vary in time 

— Time-varying heat transfer coefficients from detailed flow computations curve fit as a function of altitude, velocity, and size 


Test Case: Chelyabinsk [3] 


Heat Transfer Coefficient, C H , Model 


Basic Assumptions: 

Hyperbolic excess velocity: 

Altitude at entry: 

Relative velocity at entry: 

Relative entry angle: 

Relative heading angle: 

Geographic latitude at entry: 

Oblate rotating Earth 
Gravitational model includes J 2 term 

US-1976 atmospheric model 


15.0 km/s 

95.0 km 

19.0 km/s 
-18.5 deg 
-76.6 deg 
54.5 deg 


Meteoroid Assumptions: 

Shape: 

Density of meteoric material 
Aerodynamic model: 


Sphere 
3300 kg/m 3 
Sphere 


Sensitivity study to entry mass, heat transfer 
coefficient, heat of ablation, and fragmentation 
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Curve fit expressions are to be used for z > 15. 5 km 

C H for different velocities and diameters obtained through linear interpolation 
An example "quality of fit" plot generated with curve fit . 


Ch vs. Altitude 

V = 16 km/sec, Base Diameter = 10 m 
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Basic Plots for Variable C H and Double Fragmentation (Case 2_D) 


Trajectory Simulation Process with Meteor Physics Equations 


Entry State (t = 0) 
(Velocity, Entry Altitude, Entry 
Angles, Geographic Location, 
Meteor Shape, Size and Mass) 


Trajectory 

(Single or Multiple Bodies 
traversing over an Oblate Earth 
with J 2 harmonic and tabulated 
atmospheric model) 


Meteor Shape(s) 
(Tabulated or Newtonian 
Aerodynamic Model, Surface 
Recession, Fracture and Breakup) 


Heat Flux Environment 
(Convective and Radiative Heat 
Flux Modeled as Heat Transfer 
Coefficient, C H ) 


Mass Loss 

(Ablation, Sublimation and 
Melting) 


Mass Change Equation 
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Cross sectional area at entry 

Mass at time t 

Relative velocity at time t 

Atmospheric density time t 

Heat of ablation 

Heat transfer coefficient 


Shape Change Equation 
(assuming spherical meteor) 
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Sensitivity to Basic Assumptions: Entry Mass, Fragmentation, 

Heat Transfer Coefficient and Heat-of-Ablation 


Altitude versus Velocity Chelyabinsk Single Body Meteor 
Compared to Heat Transfer Coefficient Data Points 


Altitude versus Velocity for Chelyabinsk Fragmented Meteor 
Compared to Heat Transfer Coefficient Data Points 
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Chelyabinsk Observed 
Case 2_A 
Case 2_B 
Case 2_C 
Case 2_D 
Heat Transfer Coef. 


1 1 1 1 1 1 1 1 1 ■ 1 1 1 1 1 1 1 1 1 1 1 


□ 


□ 

□ 

□ 

□ 

□ 

□ 


□ 

□ 

□ 

□ 

□ 

□ 

□*. 


; 

i 


□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 


0 


20 


5 10 15 

Velocity, km/sec 

Case 1__A: Single body with C H = 0.1 
Case 1_B: Single body with time varying C H 

Entry mass: 13 xlO 6 kg 

Meteor radius at entry: 9.8 m 

Heat-of-Ablation: 8 MJ/kg 

Case 1_C: Time varying C H for a single body. Entry mass and heat of ablation tuned 
to overlay Chelyabinsk observations 
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Case 2_A: Single fragmentation event at 40 km, and C H = 0.1 
Case 2_B: Two fragmentation events at 40 & 30 km, and C H = 0.1 
Case 2_C: Single fragmentation event at 40 km, and C H time varying 
Case 2_D: Two fragmentation events at 40 and 30 km, and C H time varying 

For Cases 2_A - 2_D: Meteor radius: 3.5 m at 40 km alt. Revised mass: 5.93 xlO 5 kg 

For Cases 2_B & 2_D: Meteor radius: 0.7 m at 30 km alt. Revised mass: 4.74 xlO 3 kg 


Entry mass: 6.8 xlO 4 kg 

Meteor radius at entry 17m * Myriad ways to fit observations by choice of model parameters and fragmentation events 

Heat-of-Ablation: 0.85 MJ/kg 


Problem compounded by the fact that exo-atmospheric dynamical mass not known precisely 


Conclusions, Future Work and References 


Heat Transfer Coefficient vs. Altitude 


Meteor Radius vs. Altitude 
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Stagnation Point Radiative Heat Flux versus Altitude 
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Mass vs. Altitude 
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Mass vs. Altitude 
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For Case 2_D simulation fragmentation at 40 & 
30 km altitudes assumed to occur instantly 
Fragment masses tuned to overlay simulated 
trajectory on Chelyabinsk observations. 

On a scale of 40 to 90 km altitude, mass vs 
altitude trace appears to be a straight line over 
the entire mass range 

Trace is actually parabolic when mass scale is 
expanded 

Influence of C H model is insignificant if large 
changes occur in meteor mass due to 
fragmentation. 


Traj , an established trajectory simulation tool successfully modified for meteor entries 
Improvements include: 

• Simple mass loss equation of meteor physics 

• Time-varying heat transfer coefficient based on detailed flow computations 

• Ability to specify fragmentation events 

Updated version of Traj tested against Chelyabinsk observations 

Traj can now be used to establish sensitivity of trajectories to various meteor parameters 
Leaves open the issue of verification/validation of Traj and additional test cases are needed 
Could tektites [4] be used as additional test cases? 

Advantages of simulating tektite entries into Earth's atmosphere 

• Exo-atmospheric shapes are definitely spherical 

• Small sizes and (sub)orbital entry velocities 

• Problem is dominated by convective heating and melting 

• Melted shapes are aerodynamically stable 

• Chemical composition of australite tektites is statistically well defined 

• Serve as a good foundation for the tougher meteor entry problem 



Recovered australite tektite 
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